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ACOUSTIC  WAVES  INCIDENT  ON  A  SEAWATER/SEA  ICE  INTERFACE 

by 

Kenneth  C.  Jezek 

INTRODUCTION 

Acoustic  signals  propagated  in  the  Arctic  Ocean  can  scatter  at  the  ocean 
bottom  and  at  the  sea  Ice  canopy.  Crary  (1954)  studied  the  problem  of  com- 
presslonal  (P)  and  distortlonal  (shear  or  S)  waves  incident  on  a  sea  ice/ 
seawater  Interface  In  conjunction  with  seismic  studies  on  an  Arctic  ice 
island.  More  recently,  several  investigators  have  completed  straightforward 
analyses  of  the  reflectance  of  acoustic  waves  at  a  perfectly  smooth  seawater/ 
sea  ice  interface  because  understanding  sound  propagation  in  the  Arctic  Ocean 
is  Important  to  underwater  surveillance  and  detection  schemes  (Mayer  1974, 
Mayer  et  al.  1975).  In  this  paper,  I  extend  the  results  of  the  earlier  work 
to  examine  both  the  reflected  and  transmitted  waves  at  a  seawater/sea  ice  in¬ 
terface.  This  work  is  motivated  by  acoustic  experiments  in  the  sea  ice  test 
facility  at  CRREL  where  the  back-scattered  echo  strength  from  the  underside 
of  sea  ice  and  transmission  through  the  sea  ice  in  the  forward  direction  are 
being  measured. 

THEORY 

The  theory  for  computing  the  amplitude  of  reflected  and  refracted  plane 
waves  at  a  perfectly  smooth  boundary  between  a  liquid  and  a  solid  is  well 
known  (Knott  1899).  The  solution  starts  by  defining  two  potentials,  4>  and 
that  are  related  to  the  compressional  and  distortional  wave  motions.  Then, 

wave  equations  in  4>  and  ^  are  solved  under  appropriate  boundary  conditions  to 

yield  ratios  of  the  incident  and  scattered  wave  potentials. 

The  application  of  this  approach  to  a  slab  of  sea  ice  overlying  a  semi¬ 
infinite  ocean  basin  is  limited  by  several  assumptions.  First,  multiple  re¬ 

flections  in  the  slab  are  ignored  by  postulating  that  the  incident  energy  is 
a  short  pulse  for  which  the  material  properties  are  non-dispersive.  (This 
limitation  can  be  relaxed  if  interference  between  the  multiply  scattered 
waves  is  important.)  Second,  Poisson's  constant  (a)  is  taken  to  be  equal  to 
0.25.  This  minimizes  the  algebra  involved  in  computing  the  reflection  and 
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transmission  coefficients  without  seriously  distorting  the  results.  (Pois¬ 
son's  constant  0  =  0.25  implies  that  the  compresslonal  wave  velocity  is  equal 
to  the  cube  root  of  three  times  the  shear  wave  velocity.  Roethlisberger 
[1972]  has  tabulated  values  of  the  elastic  wave  velocities  in  sea  ice  that 
show  this  assumption  is  approximately  true.)  Third,  plane  wave  theory  does 
not  accurately  predict  the  amplitude  of  lateral  waves  and,  in  this  analysis, 

I  ignore  Rayleigh  waves  and  flexural  waves.  Therefore,  the  predicted  ampli¬ 
tudes  are  likely  to  be  most  in  error  near  critical  angles  (Press  and  Ewing 
1951,  Oliver  et  al.  1954,  Mayer  et  al.  1975).  Fourth,  bulk  material  proper¬ 
ties  are  taken  as  appropriate  to  the  entire  slab  of  sea  ice.  This  assumption 
is  especially  weak  near  the  ice/ water  Interface  where  there  is  the  highest 
density  of  channels  filled  with  liquid  brine.  However,  Mayer's  (1974)  re¬ 
analysis  of  Langleben's  (1970)  data  on  the  reflectivity  of  sea  ice  suggests 
that  ice  near  the  ice/water  boundary  can  support  shear  waves.  So,  I  suspect 
the  analysis  in  this  paper  to  be  valid  as  long  as  the  wavelengths  are  much 
greater  than  the  dimensions  of  brine  drainage  channels  and  the  individual 
crystals  in  which  the  channels  are  aligned. 

Solution  of  the  wave  equations  subject  to  the  above  restrictions  yields 
three  ratios  for  the  reflected  and  refracted  P-wave  and  refracted  S-wave  po¬ 
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where  p  is  the  density  of  seawater,  p'  is  the  density  of  sea  ice,  8  is  the 
shear  wave  velocity  in  sea  ice  and  c  is  the  horizontal  phase  velocity.  (c 
becomes  infinite  at  normal  Incidence.)  The  remaining  variables  are  defined 
by 
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where  0p,  8p'  and  0S  are  the  deflection  angles  of  the  reflected  P-wave, 
the  refracted  P-wave  and  the  refracted  S-wave,  a  is  the  compresslonal  wave 
velocity  in  sea  water  and  a'  is  the  compresslonal  wave  velocity  in  sea  ice 
(Ewing  et  al.  1957).  Note  that  these  angles  become  imaginary  if  the  critical 
angle  is  exceeded.  These  ratios  can  be  used  to  evaluate  wave  energies 
through  the  equations 


r,T  _  p*  r* 


Es  p*  s  ^ 
Et  '  0  r*2 


where  Er  is  the  energy  of  the  reflected  P-wave,  Ef  is  the  energy  of  the 
transmitted  P-wave,  Eg  is  the  energy  of  the  transmitted  shear  wave  and  E^ 
is  the  energy  of  the  incident  wave. 

ANALYSIS  AND  DISCUSSION 

I  used  two  sets  of  velocities  and  densities  to  compute  reflected  and 
transmitted  energies.  In  the  first  case  (Fig.  1)  I  took  the  "standard” 
values  quoted  by  Mayer  et  al.  (1975).  Specifically,  p’  =  0.917  g/cm  ,  p  = 
1.03  g/cm3,  a  3  1500  m/s,  a'  3  3500  m/s  and  8  =  1550  m/s.  The  second  set  of 
material  constants  (Fig.  2)  corresponds  to  values  l  believe  appropriate  to 
the  conditions  In  the  CRREL  sea  ice  test  facility.  These  are:  p'  =  0.91 
g/cm3;  p  3  1.03  g/cm3;  a  3  1450  m/s;  «'  3  3500  m/s;  6  3  1550  m/s.  (The  wave 
speed  in  water  was  measured  by  carefully  timing  signals  reflected  off  the 
bottom  of  the  sea  Ice  and  detected  by  a  receiver  at  a  known  depth.  The  dens¬ 
ity  of  the  sea  ice  was  estimated  by  determining  the  volume  and  weight  of 
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c.  Transmitted  S-wave.  d.  Relative  energy,  curves  for  the 

reflected  compresslonal  (P),  trans¬ 
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transmitted  shear  (S)  waves  are 
shown. 


Figure  1.  Reflection  and  transmission  coefficients,  standard  values  of  mate¬ 
rial  properties  (p'  *  0.917  g/cm3,  p  *  1.03  g/cm3,  a  =  1500  m/s,  a'  =  3500 
m/s  and  3  *  1550  m/s). 


slabs  removed  from  the  sea  Ice  sheet.*  The  compressional  wave  velocity  was 
crudely  estimated  on  the  basis  of  the  strength  of  normally  reflected  acoustic 
waves  from  the  base  of  the  sea  ice.)  Mayer  et  al.  (1975)  give  a  more  com¬ 
plete  discussion  of  the  change  in  reflection  coefficient  with  changing  mate¬ 
rial  properties. 

The  results  (Fig.  1  and  2)  are  presented  in  terms  of  the  ratios  of  the 
scattered  to  Incident  waves  expressed  in  dB.  Also  presented  are  curves  of 
the  relative  energy  of  each  wave.  Phase  information  is  suppressed  In  this 
presentation. 

*  Personal  communication  with  S.  Arcone,  CRREL,  1984. 
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Transmitted  S-wave. 


d.  Relative  energy. 


Figure  2.  Reflection  and  transmission  coefficients,  test  facility  properties 
(o'  =  0.91  g/cm3,  p  =  1.03  g/cm3,  a  =  1450  m/s,  a'  ~  3500  m/s  and  B  =  1550 
m/s). 

The  results  for  the  two  cases  that  I  considered  show  transmitted  P-wave 
cutoff  angles  at  about  25°  and  S-wave  cutoff  angles  at  about  70°.  Reflected 
P-wave  amplitudes  are  nearly  constant  to  within  a  few  degrees  of  the  cutoff 
angle.  For  Incident  angles  between  30°  and  b0° ,  the  reflected  P-wave  is 
severely  attenuated.  S-wave  amplitudes  increase  in  this  region  and  there  is 
nearly  complete  transformation  of  the  Incident  acoustic  energy  into  shear 
wave  energy  in  the  sea  ice. 

The  observation  that  a  significant  partition  of  energy  from  incident 
acoustic  waves  to  transmitted  shear  waves  occurs  at  Incident  angles  between 
30°  and  60°  suggests  that  mode  conversion  is  an  Important  factor  when  consid¬ 
ering  sound  propagation  beneath  the  Arctic  sea  ice  cover.  This  also 
suggests  that  sensors  placed  on  the  surface  of  the  sea  ice  to  detect  under¬ 
water  sound  will  be  more  likely  to  detect  P-waves  at  incident  angles  less 
than  30°  and  more  likely  to  detect  S-waves  at  angles  between  30°  and  60°. 
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